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ABSTRACT 


ERTS-1 MSS digital data from two consecutive passes over 
the Los Angeles area were analyzed in a variety of ways to determine 
the effect of haze on the data. One day, 21 Oct 1972, was quite hazy 
over Los Angeles. There was light haze over some of the test site 
and some areas were quite clear. The haze (smog) over Los Angeles 
was 3, 500 feet thick. The visibility was about 5 miles. The particu- 
late scattering optical depth was found to be 0. 30 at 0. 70 ]im from 
analysis of solar aureole brightness measurements. The second day, 

8 Nov 1 972, was quite clear with about 10% of the amount of haze present 
on 21 Oct. The data were analyzed in a variety of ways. 

The basic effect of the haze on the MSS data was found to be 
an overall increase in apparent brightness of the scene. The magnitude 
of this brightness increase was found to be approximately independent 
of surface reflectivity, over the range of reflectivities encountered in the 
test area. The magnitude of the brightness increase was found to decrease 
with increasing wavelength. The amount of haze present caused no signif- 
icant decrease in the power spectral density, even at the highest spatial 
frequencies in the MSS data. It is estimated that the total amount of 
haze can be measured over water by analysis of the apparent brightness 
in MSS bands 4, 5, and 6. Only the total amount, not the vertical distri- 
bution can be determined. The accuracy is estimated to be 1/2 to 1/3 
of the amount of haze present over Los Angeles on 21 Oct 1972. The 
amount of haze can also be determined over land with about the same 
or slightly better accuracy, by noting the variation in apparent brightness, 
if some form of ground truth for that area is known at times of previous 
ERTS passes over the area. This investigation was based on limited 
data over a single test site, however, the results are thought to be 
generally applicable. 
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INTRODUCTION 


The purpose of this investigation was to demonstrate the feasi- 
bility of and determine the optimum method for remote haze monitoring 
by a satellite, by analysis of ERTS-1 Multi-Spectral Scanner data. 

The approach was to obtain as much ground truth as possible, 
at the times of ERTS-1 passes over the Test Site (the Los Angeles area), 
in order to choose appropriate data for the numerical analysis and to 
correlate the amount of haze with the magnitude of its effect on ERTS-1 
MSS data. The type and magnitude of these effects pointed the way toward 
optimum methods for analysis of ERTS data to obtain the amount of haze 
present, and an estimate of the accuracy of these methods. 

In many respects the ERTS-1 MSS was ideally suited for an 
analysis such as this. The wavelengths from 0. 5 to 1.1 jam lie in a region 
where Rayleigh scattering is rapidly becoming less important, and where 
aerosols play a significant role in atmospheric scattering. The presence 
of four wavelength bands allows one to determine the optimum wavelength 
for detection of haze, and the relatively high spatial resolution allows one 
to determine the effects of haze on contrast at high spatial frequencies. 

On the other hand, the primary motivation for flying the ERTS-1 
satellite was to image the surface of the Earth. The sensor was configured 
to look straight down, where atmospheric effects are minimized. 
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GROUND TRUTH 


It was necessary to determine the amount and distribution of 
haze over the test area at the times of ERTS-1 passes. It was also nec- 
essary to establish the presence or absence of clouds, especially high 
thin clouds, which might perturb the analysis. These data were required 
in order to choose appropriate data for numerical analysis and to correlate 
the amount of haze with the magntidue of its effect on ERTS multispectral 
imagery. Several types of ground truth data were obtained. They are 
described in this section. The data are compiled in Appendix A. 

Four types of ground truth data were obtained. 

a) Visibility and relative humidity data from airports 

b) Air Pollution Control District data 

c) Solar aureole brightness 

d) Visual and photographic observations from a light 
airplane. Vertical distribution of haze. Vertical 
temperature profile. 

The first two were obtained from outside ongoing sources. The last two 
were actively gathered by the investigators specifically for this analysis. 

a) Visibility and Relative Humidity Data from Airports 

Horizontal visibility at ground level is determined hourly, on 
the hour, at controlled airports and broadcast continuously for the use of 
private and commercial pilots. 

ERTS-1 passed over our test site between 5 minutes before and 
one minute after the hour. Therefore, the time of visibility observations 
coincided nicely with ERTS-1 passes. 

Visibility reports were recorded at the time of every ERTS-1 
pass over our test area from the following airports, the locations of which 
are shown in Fig. 1. 
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1 Los Angeles International 

2 Burbank 

3 Van Nuys 

4 Ontario 

5 Orange County 

6 Long Beach 

7 Torrance 

8 Palmdale 

9 Palm Springs 

The visibility data was not expected to be uniquely correlated 
with the effects of haze on the ERTS data, because it is a measure of 
horizontal visibility and does not reflect the vertical profile of the haze. 
Also, it is a semi -subjective measurement which may vary from observer 
to observer. It does have the advantage that it is readily available from a 
large number of sites and is familiar to the technical community. 

The visibility data coupled with the measurements of the vertical 
structure of the haze layer made by the investigators (d) is thought to be 
a useful measure of the amount of haze. 

Relative humidity measurements, which were available from 
many of these airports, were collected to compliment the Air Pollution 
Control District particulate density data, as explained below. 

b) Air Pollution Control District Data 


The Air Pollution Control District, County of Los Angeles, 
records the amounts of various pollutants at ground level hourly at 12 
stations. The locations of these stations, designated by letters A-L, 
are shown on Fig. 1. 


A Los Angeles 

B Azusa 

C Burbank 

D West Los Angeles 

E Long Beach 

F Reseda 


G 

H 

I 

J 

K 

L 


Pomona 

Inglewood 

Whittier 

Newhall 

Lancaster 

Pasa'dena 
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The pollutants measured are: 


NO 

(parts 

no 2 

(parts 

°3 

(parts 

Hydrocarbons 

(parts 


per hundred million) 
per hundred million) 
per hundred million) 
per million) 


ch 4 

CO 

so 2 

Particulate 

Density 


(parts per million) 

(parts per million) 

(parts per hundred million) 

(Km units X 10) 


Km is a measure of particulate density based on the reduction 
of light reflected from filtered particulates compared with that reflected 
from clean filter paper. 

R 

Km = 10 log jq 

where R Q is the reflectance of clean Whatman #52 filter paper at 

0. 40 (am. 

3 

R is the reflectance after one m of air is filtered through 
, 2 

an area of one cm . 

It has been empirically determined (Ref. 1 ) that in the case of 
carbon particles: 

11 3 

Carbon particulate density « (7.7) x (Km) ’ |igms/m . 

Km measurements alone are not sufficient to define the optical 
degradation of the atmosphere, because visibility has been found to be 
strongly dependent on relative humidity for a given value of Km (Ref. 1 ). 
Evidently, liquid water forms about the particulate contaminants to an 
amount determined by the relative humidity, thus significantly increasing 
the scattering and absorption of radiation in the atmosphere. For this 
reason, relative humidity measurements made at the local airports were 
collected along with visibility data. 
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c) Solar Aureole Brightness 

A solar aureole monitor was designed, and two units were built. 
These instruments measured the brightness of the solar aureole at angles 
of 1.5, 3.0, and 6.0 deg from the center of the solar disc along the almu- 
cantar, i. e. , parallel to the horizon, in a spectral band extending from 
0.6 (jm to 1.0 14 m. 

The brightness of the solar aureole results from small angle 
scattering, principally from aerosols, and is therefore a very sensitive 
monitor of the haze producing aerosols along the path between the observer 
and the sun (Refs. 2, 3, 4). The broad spectral band was chosen to approxi- 
mate that of the sum of the ERTS-1 sensors, in order to maximize the 
correlation of these data with the effects of haze on ERTS-1 data, and to 
reduce the Rayleigh scattering contribution to the aureole measurements. 

These data proved to be less useful than orginally envisioned 
because only two predetermined locations could be covered during each 
ERTS-1 pass, whereas the test area covers many areas of vastly different 
haze conditions. Even a small cloud between the observer and the sun at 
the time of a pass rendered the aureole monitor useless. There have been 
very few measurements of the brightness of the solar aureole, therefore, 
there is almost no base with which to compare the aureole data. 

d) Visual and Photographic Observations from a Light Airplane; 

Vertical Distribution of Haze; Vertical Temperature Profile 

Flights were made over the test area, weather permitting, at 

the times of ERTS-1 passes in a light plane (Cessna 172 or 150). 

Visual observations, from above the haze layer, of the location, 
extent, and relative amounts of haze were recorded. Photographs were 
taken to document the visual observations and to aid in establishing relative 
amounts of haze from pass to pass. These photographs were taken using a 
Honeywell Pentax Spotmatic 35 mm camera with a 50 mm f/1.4 Super- 
Multi-Coated Takumar lens and a U.V. filter. The Kodachrome II film 
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was processed by Kodak. Photographs were taken looking both down and 
toward the horizon in a number of directions. Photographs were taken 
over a variety of locations in the test area. Time, location, altitude, 
direction, and camera settings (typically 1/500 Sec at f/4. 9) were noted 
and later written on the cardboard mount of each color slide. Examples 
of some of these photographs, comparing the haze conditions on 
21 October 1972 and 8 November 1972, are shown in Figs. 2 — 7. It was 
possible to cover a large part of the test area within 30 min of the ERTS-1 
passes in this manner. Even though this method was subjective, the 
investigators feel that it was the best available method for establishing 
the haze distribution over the test area at the time of ERTS-1 passes. 

The vertical distribution of haze was observed while climbing 
out from the airport. The top of the haze layer, if one existed, was deter- 
mined by flying at the top of the layer and noting the altimeter reading. 

The vertical temperature profile was taken by recording aircraft 
outside air temperature vs. altitude. 

Summary of Ground Truth Data Obtained at the 
Times of ERTS-1 Passes 

Table I lists the dates of ERTS-1 passes over our test area and 
the amount of ground truth obtained during each pass. 

Through November 1972, ground truth was obtained not only for 
the pass over Los Angeles, which is the heart of the test area, but also for 
the adjacent pass (one day earlier) over San Bernardino, which is on the edge 
of our test area. Data over San Bernardino is less useful than that over 
Los Angeles and the ground truth was more difficult to obtain. However, 
because of the unpredictability of the weather and of the lifetime of space- 
borne sensors, we thought it wise to cover all possible passes until a 
minimum amount of data useful for our analysis was obtained. 

The effort to collect ground truth was stopped on 18 July 1973, 
because it would not have been possible to obtain and analyze digital tapes 
from passes at later dates within the time frame of this contract. 

The ground truth data are contained in Appendix A. 
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Table I. Summary of Ground Truth Data Obtained at the Time of ERTS-1 Passes 


Date 

Area of 
ERTS-1 Data-’' 

Data 

Obtained** 

Notes 

9 Aug 72 

10 Aug 72 

SB 

LA 

a) b) 
a) b) 

We were not aware of these 
ERTS-1 passes, so no active 
ground truth data was obtained. 

27 Aug 72 

SB 

b) 

Very cloudy 

28 Aug 72 

LA 

b) 

Too cloudy to obtain active 
ground truth 

14 Sept 72 

SB 

a)b)c)d) 

Cloud free. Heavy haze over 
part of the area, other parts 
were clear 

15 Sept 72 

LA 

a)b) 

Too cloudy to obtain active 
ground truth 

2 Oct 72 

SB 

a)b)d} 

Too cloudy for analysis 

3 Oct 72 

LA 

a)b) 

Very cloudy 

20 Oct 72 

SB 

a)b) 

Very cloudy 

21 Oct 72 

LA 

a)b)c)d) 

Essentially cloud free. Heavy 
haze over part of the area, other 
parts were clear 

7 Nov 72 

SB 

a)b) 

Very cloudy 

8 Nov 72 

LA 

a)b)c)d) 

Essentially cloud free. Clear 
day. Ideal for comparison with 
21 Oct 72 data 

25 Nov 72 

SB 

a)b)d) 

High cirrus clouds, very little 
haze 

26 Nov 72 

LA 

a)b)c)d) 

Clear day 

14 Dec 72 

LA 

a)b) 

Cirrus clouds over test area 

1 Jan 73 

LA 

a)b) 

Clear with some dust in the air. 
Holiday 

19 Jan 73 

LA 

a)b) 

Partly cloudy, very little haze 

6 Feb 73 

LA 

a)b) 

Very cloudy 

24 Feb 73 

LA 

a)b) 

Very cloudy 

14 Mar 73 

LA 

a)b)c)d) 

Very clear day 

1 Apr 73 

la 

a)b)c)d) 

Partly cloudy, very light haze 

19 Apr 73 

LA 

a)b)c)d) 

Scattered clouds, light to 


moderate haze 



Table I (Cont. ) 


7 May 73 

LA 

a)b)c)d) 

Scattered low clouds , moderate 
to heavy haze, some high cirrus 
clouds 

2 5 May 73 

LA 

a)b) 

Very cloudy 

12 June 73 

LA 

a)b) 

Very cloudy 

30 June 73 

LA 

a)b) 

Very cloudy 

18 Jul 73 

LA 

a)b)d) 

Broken clouds, moderate to 
heavy haze, high cirrus clouds 


ERTS-1 pass over San Bernardino SB 

Los Angeles LA 

a ) Visibility and relative humidity data from airports 

b) Air Pollution Control District data 

c) Solar aureole brightness 

d) Visual and photographic data 
Vertical distribution of Haze 
Vertical temperature profile 
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digital processing techniques 


AH quantitative analysis of the ERTS data has been done from 
Computer Compatible Tape (CCT) sets provided by NASA. The basic 
program for reading these tape sets at Aerospace was made as flexible 
as possible. The program is modular, separating into various subroutines 
the different operations on the tapes, so that changing analysis procedures 
would require as little recoding as possible and therefore be as simple 
as possible. All of these programs were written in Fortran for the 
Control Data Computers at The Aerospace Corporation. Because they 
employ properties of the Aerospace operating system which are probably 
not generally available the programs themselves are probably not of 
general interest, although we would be glad to supply them, as is, to 
anyone who is interested. Because of this uniqueness to the Aerospace 
Corporation computer environment we describe what the programs do 
rather than details of how the programs are controlled and used. We 
discuss first the basic program and then the various analysis procedures 
developed for use with the basic program. 

Basic Computer Program 

The flow of operations in the basic computer program is 
shown in Figure 8. The names of the subroutines performing the various 
operations are indicated under the lower right corner of the box containing 
the operations description. 

The spatial region to be analyzed is selected by cards read 
by the routine DESLOC. This is specified by giving the beginning and ending 
scan line numbers in the frame and the beginning and ending point numbers 
in the scan line. This routine also recognizes the stop directive on the 
input cards. 

The analysis part of the program consists of 3 subroutines 
(or 3 entry points into one subroutine). The names of the entry points are 
ANPAR, LINAN and FINAN. The entry point ANPAR is called before the 
analysis begins and permits the input of various parameters needed for 
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the analysis. The program then begins to deliver, to the LINAN entry point 
one scan line at a time, the MSS data for all four bands. After all the 
desired lines have been delivered and analyzed, the entry point FINAN is 
called to do any summarizing needed to complete the analysis procedure. 

DESLOC is called again after the call to FINAN. An attempt 
is made to input another spatial location and the whole procedure is repeated. 
Since the data is on the tapes in sequential scan line order it is much more 
efficient to order the locations desired for analysis in ascending order of 
line number. The program will handle out of order requests correctly 
but at the expense of rewinding the input tape files and skipping records 
which would not otherwise be necessary. 

The difference in the computer word format between the XDS 
computers at the ERTS NDPF and the CDC computers available at Aero- 
space made the writing of these programs somewhat more difficult than 
it might have been otherwise. Since the XDS format is essentially the 
IBM 360-370 format it was deemed desirable to create a somewhat 
more general purpose routine for the conversion of data from XDS format 
to CDC format, than might have been required strictly for the ERTS 
analysis. This routine is described in Appendix C. 
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RESULTS 


On the basis of ground truth information, the data from the 
ERTS-1 passes on 21 Oct 1972 and 8 Nov 1972 were picked for extensive 
digital analysis. The ’’observation identifier” for these frames are 
1090-18012 and 1108-18014. During the 21 Cct 1972 pass much of the 
test area was quite hazy. (See Fig. 2, 4, 6 and Appendix A). There 
were also regions of moderate and light haze as well as clear areas. 

At the time of the 8 Nov 1972 pass, most sections of the test area, in- 
cluding regions that were hazy on 21 Oct, were quite clear (see Fig. 3, 

5, 7 and Appendix A). We had complete ground truth for both passes, and 
both were essentially cloud free. Differences in terrain features and 
solar elevation angle are minimal, because these were consecutive passes. 
These images were compared in various ways to determine the effects of 
haze on ERTS-1 images. 

Quantitative Assessment of Haze 

A quantitative assessment of the amount of haze present at 
various locations on 21 Oct 72 and 8 Nov 72 from the ground truth data 
listed in Table A- II, pp A- 13 and A-15. 

First, the particulate density measurements. Km, made at 
ground level by the Air Pollution Control District were studied. In an 
attempt to determine the relative importance of solid particulates and 
water droplets, all available data from Table A- 11 were used to make a 
least squares fit to the form 

visibility (miles) * 

to determine the parameters a and b, R.H. is the relative humidity, in 
percent. The data from Table A -II, cover ranges of 29 to 88% relative 
humidity, 2. 5 to 40 miles in visibility, and Km values from 0. 9 to 4. 5. 

The best fit values were 


a = 927 miles 
b = 0. 10 
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b is negligibly small compared with reasonable values of relative 
humidity. Thus the atmospheric scattering becomes very small when 
either Km or the relative humidity becomes very small. This implies 
that water droplets are the primary scattering material and that solid 
particulates play the role of condensation centers. We don't see any 
way to determine aerosol concentration from the Km measurements. 

The solar aureole brightness measurements may be used to 
obtain aerosol concentrations; however, some assumptions about 
materials and size distributions must be made . For reasons discussed 
in the preceding paragraph, water was chosen as the material. 

Deirmendjian (Ref, 5) has calculated the scattering of light of wave- 
length 0, 70^m (near the peak response of the aureole monitor, 0,75 ^m) 
by water droplets. He made these calculations for two particle size 
distributions M (marine) and L (continental). The normalized total 
phase functions are shown as a function of scattering angle in Fig. 9. 

The relative aureole brightness values at 1 1/2, 3, and 6 deg, measured 
on 21 Oct 72 agree reasonably well with the L-type scattering phase func- 
tion. Therefore the following calculations will be made using Deirmendjian' s 
L-type particle size distribution (Ref. 5). 

Following the nomenclature of Ward et al. (Ref. 4), 


b t (x,4») = 


where: 



e = 
4 * 

x = 
R(4) = 


P<4,x) = 
t r (\) = 
T p (\) = 

t(\) = 


H(X) exp[-T(X)sece][T R (X)R(40 + r p (X)P(4, X)] sec 0 


Aureole brightness 

unattenuated solar brightness (same units as B,, e. g. , 

Watts /cm^. S r • (j,m) t 

solar angle from the zenith 

scattering angle 

wavelength 

the normalized Rayleigh scattering phase function = 

3/8 (1 + cos^ iJj) 

The normalized particulate scattering phase function (Fig. 9) 
The Rayleigh optical depth, . 0292 at 0.75 (j,m (Ref. 6) 

The Particulate scattering optical depth. 

Total optical depth. 
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The values of obtained in this way were: 

r (21 Oct 72) = 0. 30 

P 

r (8 Nov 72) = 0.03 
P 

From these values of particulate optical depth, the total 

2 

number of particulates in an atmospheric column of one cm area, 

can be obtained from Deirmendjian 1 s calculations of the extinction 
coefficient (Ref. 5) (still with L-type particle size distribution). The 
results are: 


N^(21 Oct 72) = 7. 6 x 10^ particles/cm^ 

N^(8 Nov 72) = 0. 8 x 10^ particles/cm^ 

These particle number data can be used to calculate the expected 
visibility. The investigators observed that the haze on 2 1 Oct 72 was 
composed of a layer of approximately uniform density and about 3, 500 
feet thick. If the total amount of haze on 21 Oct is spread uniformly along 
a 3, 500 ft. thick layer, the number density is 

n(21 Oct 72) = 710 particles/cm^. 

The extinction coefficient p can again be obtained from Deirmendjian* s cal- 
culations, now for visible light (0.45 ^m), instead of the extinction coef- 
ficient calculated at 0. 7 p,m used to obtain number density from the aureole 
data. 


P = 0. 34 km" 1 

From this value of the extinction coefficient, ft, the visibility can be cal- 
culated (Ref. 7) as, 

visibility = 3 * 912 
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Thus, calculated visibility for 21 Oct 72 is found to be 11*5 km = 7 miles* 
This is thought to be quite good agreement with the reported visibilities 
of 5 and 6 miles. 

The ERTS-1 Multi -Spectral Scanner Digital Data 

The key characteristics of the digital data are listed in 
Table II. The analysis was carried out in units of the lease significant 
bit of the digital data word. They may be converted to radiometric units 
by multiplication with the values given in the last column of Table II. 

Images 

The first analysis package allowed the investigators to pro- 
duce images of selected areas of a frame for a particular MSS band. The 
images were produced by the Calcomp 83 5 film plotter associated with 
the CDC 7600 computer. Thirty-two gray shades were available; 35mm 
Panatomic-x film was used. Four minutes in D-96 at 70 °F was found to 
yield optimum processing. These images were somewhat skewed because 
the effect of the rotation of the earth, as the MSS scanned the scene, was 
not removed. Stripes were placed at certain scan lines in order to pre- 
cisely locate these scan lines geographically. 

This pictoral capability was developed in order to view the 
effect of haze on the image, to precisely geographically locate the pixels 
included in the various digital analyses, and to display the results of 
numerical manipulation, e. g., ratios or differences of two frames of 
data. 

Figures 10 and 11 show photographs of an area in the western 
part of the Los Angeles basin made by the investigators from MSS Band 5 
digital tapes. The area was quite hazy at the time of Fig. 10 and quite 
clear at the time of Fig. 11. The exposure times were slightly different, 
and a slightly higher contrast grade of paper was used to make the print 
of the hazy day, but the image quality appears quite comparable. Similar 
results were obtained for MSS bands 4, 6, and 7. 

Line Plots of Intensity 

Plots of the intensity of all four MSS bands along a large number 
of selected lines were made from the 21 Oct frame and the 8 Nov frame. 
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Table II. Description of ERTS-1 MSS Digital Data 


Channel 

Designation 

Spectral 

Bandpass 

( U m) 

MSS 4 

0. 5 - 0. 6 

MSS 5 

0.6 - 0.7 

MSS 6 

0.7 - 0.8 

MSS 7 

0.8 - 1. 1 


Number of Gray Shades 
(linear with brightness) 

128 (7 bits) 

128 (7 bits) 

128 (7 bits) 

64 (6 bits) 


Radiometric Equivalent 
of Least Significant Bit 
(watts/cm^« sr) 

1.95 x 10" 5 
1.57 x 10" 5 
1.39 X 10' 5 
7.30 x 10“ 5 
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These line plots were quite useful in precisely geographically locating 
picture elements along scan, as well as in examining the data to deter- 
mine the effects of haze. Figures 13 and 14 show examples of the plots 
for part of one line from the 21 Oct frame and the corresponding line 
from the 8 Nov data. There appears to be very little difference between 
the line plots for the two days except that 21 Oct data (hazy day) appear 
to be shifted up by a nearly constant amount. The amount of this shift 
decreases as wavelength (i. e., MSS band number increases. This shift 
seems to disappear when areas that were clear on both days are 
compared. 

I ntensity Value Distributions 

Brightness value distributions were computed for several areas 
from the 21 Oct frame and for the same geographical areas from the 
8 Nov frame. Figure 15 shows a histogram plot of the brightness values 
of one area from both the 21 Oct and 8 Nov data for MSS band 4. 

Figures 16-18 are similar plots of MSS band 5, 6, and 7 data. The 
area for which data is shown in Fig. 15 - 18 includes ocean, beach, 
harbor, vacant land, city, and parkland. It was hazy over this area on 
21 Oct and clear on 8 Nov. This area lies between the top pair of white 
lines in Figs. 10, 11, and 12 and extends from the ocean to about the 
prominent white streak on the right of these pictures formed by the Los 
Angeles river. These histograms show a double peak. The brighter 
peak is associated with terrain, the darker peak with water. 

The distribution is clearly shifted toward brighter values on 
21 Oct. The amount of this shift decreases with the increasing wave- 
length. It appears that the shift is greater, in absolute value, for the 
brighter peak than for the darker one. 

The fine structure in these histograms, i. e. , dips and humps 
in adjacent intensity values, must be discounted. This is because the 
ERTS digital data has missing values. In the range of brightness values 
say 20-45, about one half of the values are unpopulated (forbidden) in a 
certain line of data. The particular missing numbers are correlated 
with a particular detector, that is, they repeat every six lines. About 
the same number of values are missing for the other detectors, but the 
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particular values which are forbidden vary from detector to detector. 

The histograms shown in Figs. 15-18 are averaged over eleven 
adjacent lines. However, these are values which are forbidden for 
as many as five of the six detectors. These values show dips in the 
histograms. Other values are forbidden in as few as one detector; 
these values show peaks in the histogram plots. This phenomena was 
not present in MSS band 7. 

Mean values of brightness of both water and terrain were 
computed from the data shown in Figs. 15 - 18. This area includes 
the city of Inglewood (Location H on Fig. 1) and is referred to as the 
Inglewood area. Similar mean values were derived from brightness 
distributions around San Pedro (between locations 6 and 7 on Fig. 1) 
including part of the harbor area, and Palmdale (just south of 
Location 8 on Fig. 1), where it was clear on both 21 Oct and 8 Nov. 

The terrain data are tabulated in Table III. The mean values of bright- 
ness over water are shown in Table IV. The differences in the mean 
values. A, is also listed. The high apparent brightness of water in 
MSS band 4, and to a lesser extent in band 5, even during clear times 
is caused by reflection of the incident sunlight into the ERTS sensor by 
the clear atmosphere (Rayleigh scattering) rather than a high reflectivity 
of water at these wavelengths. The Rayleigh scattering into the ERTS 
sensor is nearly independent of solar elevation angle over the range of 
elevation angles encountered in this investigation. 

Table IV shows that there was a significant increase in the 
apparent brightness of water in MSS bands 4, 5, and 6 during hazy con- 
ditions at Inglewood and San Pedro, caused by the haze scattering incident 
sunlight toward the ERTS sensor. Palmdale, which was clear on both 
days, showed about the same apparent brightness for water on the two 
days. The magnitude of this effect decreases with increasing wavelength. 
For example, at San Pedro, the A's for MSS bands 4, 5, 6, and 7 are: 

5. 8, 4. 8, 3.4, and 0.4, in units of the least significant bit of the ERTS 
MSS data words. When converted to radiometric units using the data from 
Table II, these A’s become 11.3, 7.5, 4.7, and 2. 9 x 1 0 -5 watts/cm 2 , sr. 
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Table III. Mean Values of Brightness 

jjj 

Distributions Over Terrain 


Inglewood 

21 Oct 72 
(Hazy) 

8 Nov 72 
(Clear) 

A 

MSS Band 4 

35. 6 (35. 1) 

27.7 

7. 9 (6.4) 

5 

30.9 (29.7) 

24.7 

6.2 (5.0) 

6 

31.4 (30.5) 

26. 1 

5. 3 (4.4) 

7 

13.4 (13.2) 

12.0 

1.4 (1.2) 

San Pedro 

21 Oct 72 
(Hazy) 

8 Nov 72 
(Clear) 

A 

MSS Band 4 

35. 5 (33.8) 

29.0 

6.5 (4.8) 

5 

30.8 (29. 5) 

25.3 

5. 5 (4.2) 

6 

25.9 (25. 1) 

22.7 

3.2 (2.4) 

7 

9.4 (9. 3) 

8. 7 

0.7 (0.6) 

Palmdale 

21 Oct 72 
(Clear) 

8 Nov 72 
(Clear) 

A 

MSS Band 4 

39. 7 (37.0) 

37.0 

2.7 (0) 

5 

43.8 (41.5) 

41.5 

2.3 (0) 

6 

40. 6 (39.2) 

39.2 

1.4 (0) 

7 

17.7 (17.4) 

17.4 

0. 3 (0) 


Values in parentheses are the brightness values after normalizing the 
same solar elevation angle as 8 Nov 72, The observed variations in the 
Palmdale data were used to make the corrections at Inglewood and San 
Pedro. 



Table IV. Mean Values of Brightness 
Distributions Over Water 


Inglewood 

21 Oct 72 
(Hazy) 

8 Nov 72 
(Clear) 

A 

MSS Band 4 

21.0 

15.8 

5.2 

5 

10. 5 

5. 9 

4.6 

6 

5. 9 

2.5 

3.4 

7 

0.3 

0 

0. 3 

San Pedro 

21 Oct 72 
(Hazy) 

8 Nov 72 
(Clear) 

A 

MSS Band 4 

23.4 

17. 6 

5.8 

5 

12. 1 

7.3 

4. 8 

6 

6.6 

3.2 

3.4 

7 

0.4 

0 

0.4 

Palmdale 

21 Oct 72 
(Clear) 

8 Nov 72 
(Clear) 

A 

MSS Band 4 

18. 1 

17. 9 

0.2 

5 

9.3 

10,2 

-0.9 

6 

4.4 

5. 9 

-1.5 

7 

0 . 1 

0.4 

-0. 3 
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If bands 5, 6, and 7 are normalized to the same amount of incident 
sunlight as band 4, the relative values of these A’s become 11. 3, 8. 9, 

7. 1, and 2. 3. The haze, therefore, is monotonically less effective in 
scattering incident sunlight into the ERTS sensor as the center wavelength 
of the band becomes larger. This is a reasonable behavior. Deirmenjian' s 
calculations of scattering by water droplets with an L-type particle size 
distribution shows a decrease of the extinction coefficient from 0. 04797 km 
to 0, 03953 km When the wavelength is increased from 0.45 to 
0.70 |jim (Ref. 5). (These extinction coefficients are calculated for 100 
particle s /cm . ) Plass and Katlawar (Ref, 9) have calculated the amount 
of sunlight scattered away from the earth for different amounts of haze at 
wavelengths of 0. 7 and 0. 9 M,m for an ocean surface. They find less 
backscatter per unit incident sunlight at 0. 9 p,m than at 0,7 p,m, but not 
nearly as much less as we observe in comparing MSS bands 6 and 7, The 
results then look quite reasonable, except that the observed A in band 7 
appears anomalously small when compared with bands 4, 5, and 6. 

Table III shows that there was an increase in apparent brightness 
of terrain at Inglewood and San Pedro on the hazy day. These A also de- 
crease with increasing wavelength in a reasonable way, except for band 7 , 
which again shows a smaller A relative to the other bands than expected. 

A correction for solar elevation angle should be made in com- 
paring terrain brightnesses on these two days. The magnitude of this 
variation, as well as that caused by variations in the terrain, has been 
minimized by analyzing data from consecutive passes. It is not possible 
to accurately made this correction for a scene containing such a variety 
of shapes and materials; however, as the correction is small, an approxi- 
mate correction should suffice. The Palmdale area was clear on both days, 
therefore, it was possible to use the observed variation in terrain bright- 
ness at Palmdale to approximately correct the Inglewood and San Pedro 
data. The corrected apparent brightness values and the associated A r s 
are shown in parenthesis in Table III. 

We have here a very significant result. That is, haze causes 
an overall increase in scene brightness. A# which is approximately 
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independent of scene reflectivity over typical (not highly reflective) 
terrain. The magnitude of A decreases with increasing wavelength. 

Minimum Brightness Values 

One of the easiest parameters to determine from a set of 
data is the apparent brightness of the picture element with the smallest 
apparent brightness of the set. Table V shows the results of such an 
analysis of the same areas treated in the previous section. These 
values can be determined from the histograms shown in Figs. 15-18 
for the Inglewood area. The advantage of this analysis is that the 
entire histogram need not be produced or analyzed to obtain the minimum 
brightness value. 

The values of the differences, A's, of minimum brightness 
values on the hazy and the clear days, listed in Table V, are seen to be 
nearly equal to the mean brightness differences over water and terrain 
(Tables III and IV). for bands 4 and 5. 
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Table 


Inglewood 

MSS Band 4 

5 

6 
7 

San Pedro 

MSS Band 4 

5 

6 
7 

Palmdale 

MSS Band 4 

5 

6 
7 


V. Minimum Populated Values of 
the Brightness Distribution 


21 Oct 72 8 Nov 72 A 

(Hazy) (Clear) 

19 13 6 

9 4 5 

3 1 2 

0 0 0 

21 Oct 72 8 Nov 72 A 

(Hazy) (Clear) 

20 15 5 

10 5 5 

3 1 2 

0 0 0 

21 Oct 72 8 Nov 72 A 

(Clear) (Clear) 

16 16 0 

8 9_1 

3 4-1 

0 0 0 
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Power Spectral Density 

The power spectral density was computed as a function of spatial 
frequency over the Inglewood area (same geographical area used in the 
brightness distribution) from the 21 Oct and the 8 Nov data for all MSS 
bands. For each band for every line selected the data was reformated 
into the correct form for input to a fast Fourier transform program orig- 
inally developed at Aerospace for interferometric spectroscopy (Ref. 10). 
This reformated data was copied to tape and then used as input to the Fourier 
transform program which was used to compute the power spectrum. The 
program computed the power spectral density, P(v) as a function of spatial 
frequency, from the I(x), the difference between the observed intensity and 
the average intensity on that line, x is distance along the scan line. 


P(v) = 2 a x I(nAx) e ^^ TTn ^ x 

n=-N 

The units of P are (intensity units /km ^). The intensity units in this case 
are again the magnitude of the least significant bit in the ERTS-1 MSS digital 
data tapes. 

The power spectral density was computed for the appropriate 
segment of each line, then averaged over the 11 lines included in the Ingle- 
wood area. The results are plotted in Figs. 19-22. 

The results were essentially negative. The power spectral den- 
sity may be slightly higher on 8 Nov than 21 Oct for MSS band 4 (Fig. 19), 
however, for all other bands the two days look the same. We conclude 
that there is essentially no difference in power spectral density between the 
hazy 21 Oct data and the clear 8 Nov data. 

This result is consistent with an approximately constant differ- 
ence in brightness in the data for these two days, because a baseline shift 
has no effect on power spectral density except to the zero spatial frequency 
point, which is not included in our calculation. 
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The motive for .making the power spectral density calculations 
was to look for a possible effect of small angle scattering of light by haze 
after it has been reflected by the ground toward the ERTS sensor. Such 
small angle scattering causes light to appear as if reflected from a point 
on the ground which is slightly displaced from the actual point of reflection* 
This effect reduces resolution, contrast and power spectral density for 
high spatial frequencies* This effect becomes quickly less important for 
large spatial frequencies, because scattering from aerosols is strongly 
forward peaked (Fig. 9)* Our analysis shows that for the amount of haze 
present over Inglewood on 21 Oct 72, for a nadir looking sensor and for the 
maximum resolution of the ERTS Multi -Spectral Scanner this small angle 
scattering of the reflected light is negligible. This effect is, of course, 
minimized in the case of the ERTS-1 Multi- Spectral Scanner, by its nadir 
looking geometry, because sunlight reflected from terrain passes through 
a minimum amount of atmosphere before reaching the sensor. 

Contrast Reduction 

Image contrast has been defined in many different ways, as the 
difference or ratio of intensities 1^ and 1^, normalized in a variety of ways, 
between two points in a scene. For example 



r - 11 ~ l l - R - l 

4 Ij + I 2 R + l 

The power spectral density, described in the previous section, 
is an elegant way of describing the unnormalized difference contrast as a 
function of spatial frequency. 
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Photographic contrast usually involves the ratio of two intensities, 
because the logarithmic response of film results in a constant density differ- 
ence for a given ratio, R, otherwise independent of the magntidues of Ij and 

V 

is usually used in the literature concerning contrast reduction 
by the atmosphere, where a target intensity 1^ is treated against a uniform 
background The contrast reduction is found to be a function of I 2 for a 

given atmosphere (Ref. 11). Of course, the real world seldom provides a 
constant background intensity. 

We have found that the most significant effect of haze on ERTS-1 
MSS data is an overall increase in brightness. A, which is approximately 
independent of scene reflectivity (over typical, not highly reflective terrain). 
The magnitude of A was found to decrease with increasing wavelength. 
is not changed by this effect, whereas the magnitude of C^, and are 
all decreased. For example, 


i 



1 i + A " h " A 

. ! 2 + A 


r i 

I 2 + A 


The constrast C^, of mean harbor water against mean city back- 
ground is listed for the San Pedro area from both the 21 Oct and 8 Nov data 
in Table VI. 

It appears fruitless to pursue this line of approach in this investi- „ 
gation. There is reduction in contrast, C^, C^, or C^, caused by haze. 

But there is no apparent gain in monitoring one of these contrast functions 
rather than the overall increase in brightness, A, which is the primary 
cause of the contrast reduction. 

Image of the Difference of Two Frames 

The analysis of the 21 Oct and 8 Nov data has shown that the 
principal effect of haze on the ERTS-1 MSS data is an overall increase in 
brightness, which is approximately independent of terrain reflectivity. 
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Table VI. Contrast Between Mean Apparent Water Brightness 
and Mean Apparent Terrain Brightness Over San Pedro 
for 21 Oct 72 and 8 Nov 72 


MSS Band 4 


MSS Band 5 


MSS Band 6 


MSS Band 7 


21 Oct 72 
(hazy) 


8 Nov 72 
(clear) 

23.4 - 35. 5 

-0. 341 

17.6 - 29.0 

35.5 

29.0 

12. 1-30.8 

■0. 607 

7.3 - 25.3 

30.8 

25. 3 

6.6 - 25.9 

•0. 745 

3.2 - 22.7 

25.9 

22. 7 

0.4 - 9.4 

•0. 957 

0.0 - 8. 7 

9.4 

8. 7 


o 


-0. 393 


-0.711 


-0.859 


- 1.00 


The mean apparent brightness intensities are from Tables 

III and IV. 
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This suggests a method of accurately locating the geographical extent of 
haze by imaging the difference fo two frames. 

The data of 8 Nov was subtracted from the 21 Oct data on a point 
by point basis, and an image of the difference was made in such a way that 
zero difference was imaged as neutral gray. We expect this image to show 
the geographical extent of the haze on 21 Oct as slightly lighter feature. 

Fig. 23 shows the image of the difference for MSS Band 4. The light region 
® ver land corresponds quite well with the geographical extent of haze on 
21 Oct as determined by the investigators. There is no ground truth data 
concerning the haze distribution over the ocean to compare with the sub- 
tracted image. An area of about 45 x 100 miles is covered in Fig. 23. 

The most prominent features are clouds. Those present on 21 Oct are 
white, those on 8 Nov are black. The coastline is visible because the surf 
was heavier, i. e. , brighter, on 8 Nov than 21 Oct. 

The grid of picture elements on one pass will generally not 
exactly geographically overlap the grid obtained on another pass. There- 
fore, the data from 21 Oct and 8 Nov were geographically registered to the 
nearest line number and element number at several points. Then the data 
were smoothed by averaging each 10 x 10 array of picture elements. This 
produced a lower resolution image, but one that was registered much more 
accurately relative to picture element size. The smoothed data were then 
subtracted on a point by point basis. 
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SUMMARY AND CONCLUSIONS 


ERTS-1 MSS digital data from two consecutive passes over the 

Los Angeles area were analyzed in a variety of ways to determine the effect 

of haze on the data. One day, 21 Oct 1972 was quite hazy over Los Angeles, 

but clear in other parts of the test site. The haze (smog) over Los Angeles 

was 3, 500 feet thick. The visibility was about 5 miles (the visibility must be 

3 miles or greater in order for aircraft to operate under visual flight rules). 

The particulate scattering optical depth was found to be 0. 30 at 0.7 |jm, 

from analysis of solar aureole brightness measurements. This corresponds 

5 2 

to a number of 7.6 x 10 particles/cm in a total column of atmosphere, 
assuming Deirmendjian' s (Ref. 5) L-type particle size distributions. The 
amount of haze on the clear day 8 Nov 1972, was about 1/10 of the amount 
present over Los Angeles on 21 Nov. 

The most pronounced effect, indeed the only independent effect, 
of haze on the ERTS MSS data was found to be an overall increase in appar- 
ent brightness of the scene, caused by scattering of incident sunlight toward 
the ERTS sensor by the haze* The magnitude of this brightness increase 
was found to be approximately independent of surface reflectivity, over the 
range of reflectivities encountered in our test area. The magnitude of this 
increase in scene brightness was found to decrease with increasing wave- 
length. The absolute values of this brightness increase were 11.3, 7. 5, 

-4 2 

4.7, and 1.0 x 10 watts/cm • sr*|jm for MSS bands 4, 5, 6, and 7, 
respectively. When normalized to equal amounts of incident sunlight, the 
relative values become 11.3,8.9, 7.1 and 2.3. 

We found no decrease in power spectral density, even at the 
highest spatial frequencies in ERTS MSS data, caused by the amount of 
haze present on 21 Oct 1972. 

The optimum method for analysis depends upon one's require- 
ments. We next discuss two methods, which seem feasible as a result of 
this investigation. 
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One can determine the amount of haze, using only the data from 
a particular pass, over water. The scene reflectivity is quite low. The 
Rayleigh scattering contribution is constant. Thus one can determine the 
total amount of haze from the apparent brightness in bands 4, 5, and 6. 
Many more measurements correlating measured amount of haze with 
apparent water brightness, and calculations such as those of Plass and 
Kattawar (Ref. 9) at the ERTS wavelengths, would be necessary before 
the technique could become operational, but it is nearly certain that it 
could be done. The question is, how accurate could this method be? It 
is not possible to make a statistical analysis from such a limited number 
of observations, but some estimates can be made. The accuracy will 
probably be limited by systematic errors, rather than detector noise or 
sensor calibration, etc. The systematic errors would include variations 
in the reflectivity of water, caused by sediments, waves, foam, etc. , 
and by perhaps such parameters as the size of the body of water and the 
reflectivity of the surrounding terrain. From Table IV we see that the 
ratio of the variation of apparent water brightness, for different bodies 
of water, on the haze-free day to the increase in brightness caused by 
the haze of 21 Oct is 0.38, 0. 91, 1.00 for bands 4, 5, and 6, respectively. 
Thus we conclude that the absolute accuracy of the total column of haze 
would be about 1/2 or 1/3 of the amount of haze present on 21 Oct 1972. 

If one is willing to follow the brightness history of an area 
for some time, and if he has some kind of retroactive ground truth to 
establish brightness distributions vs. amount of haze, then a near real 
time assessment of the amount of haze would be possible from analysis 
of the brightness distributions of the area in bands 4, 5, and 6. Our 
data show more consistency in A brightness caused by haze, than in 
absolute brightness over the water. It can also be applied to areas 
where there is no water, because the increase in scene brightness was 
found to be independent of scene reflectivity over a reasonably wide 
range of reflectivity. Of course one would concentrate on the variation 
of the areas of lowest reflectivity because variations caused by changing 
solar elevation angle, etc. , would be minimized. A similar technique 
is the determination of minimum brightness value for a given area. The 
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variation of the minimum brightness was found to agree quite well with 
the shift of mean brightness values for bands 4 and 5* The accuracy 
would be degraded somewhat, by looking only at minimum brightness, 
rather than the entire distribution, but the required processing time 
would be significantly reduced. 

The effects of haze would be more pronounced, and therefore 
the sensitivity of both of the above techniques increased, if the sensor 
looked several degrees away from the nadir, rather than straight down 
as it does. 

This analysis was based on limited data. However, the results 
are thought to be generally applicable. 


- 30 - 



REFERENCES 


S. R. Hall, Anal. Chem. 24, 996 (1952). 
Reiner Eiden, Tellus, XX, 380 (1968). 


A.E.S. Green, A. Deepak, and B. J. Lipofsky, Appl. Opt. 10_, 
1263 (1971). 

Gray Ward, K. M. Cushing, R. D. McPeters, A. E. S. Green, 
Appl. Opt. 12 , 2585 (1973). 

D. Deirmendjian, Electromagnetic Scattering on Spherical 
Polydispersions (American Elsevier Publishing Company, 1969), 
pp. 75, 78, 152, 167. 

Shea L. Valley, ed. , Handbook of Geophysics and Space 
Environments (Air Force Cambridge Research Laboratories, 
1965), p. 7-27. 

Seibert Q. Duntley, J. Opt. Soc. Am. 38 , 179 (1948). 

ERTS Data Users Handbook: NASA Document 71SD4249, 

Page G-14 (revised September 1972). 

Gilbert N. Plass and George W. Kattawar, Appl. Opt. 1 1 , 

1598 (1972). 

C. M. Randall, "Users Guide to Computer Programs for 
Fourier Spectroscopy," The Aerospace Corporation Technical 
Report TR -0073(9260 - 0 1) -8, 30 November 1972. 

Seibert Q. Duntley, Appl, Opt. 3_, 551 (1964). 

Max Born and Emil Wolf, Principles of Optics (Pergamon 
Press, 1959), pp. 394-397. 



FIGURE CAPTIONS 


Fig. I 
Fig. 2 
Fig. 3 
Fig. 4 
Fig. 5 
Fig. 6 


Fig. 7 


Fig. 8 
Fig. 9 


Map of the Test Area showing airports where visibility and 
humidity data were obtained (1-9), and the locations where 
air pollutants were monitored (A-L). 

The Virginia Country Club, Long Beach (near location 6 on 
Fig. I). Looking NNE from 4, 500' altitude, 43 minutes before 
the pass of ERTS-1 on 21 Oct 72. 

The Virginia Country Club, Long Beach (near location 6 on 
Fig. 1). Looking ENE from 4,500' altitude, 30 minutes before 
the pass on ERTS-1 on 2 1 Oct 72. 

Hollywood Park, Inglewood (between locations 1 and H on 
Fig. 1). Looking ENE from 4, 500' altitude, 30 minutes before 
the pass of ERTS-1 on 21 Oct 72. 

Hollywood Park, Inglewood (between locations 1 and H on 
Fig. 1). Looking NE from 4, 500' altitude, 9 minutes before 
the pass of ERTS-1 on 8 Nov 72. 

Looking ENE from an altitude of 4, 500' over Santa Monica 
(near location D on Fig. 1) 25 minutes before the ERTS-1 
pass on 21 Oct 72. The distance to the nearest mountain 
range is 23 miles. 

Looking NE from an altitude of 3, 500' over Culver City 
(between locations 1 and D on Fig. 1) 5 minutes before the 
ERTS-1 pass on 8 Nov 72. The distance to the nearest 
mountain range, the same as that shown in Fig. 8, is 26 
miles . 

I 

Flow chart of the ERTS-1 program. 

Normalized, total phase functions for scattering from water 
droplet haze with both marine and continental particle size 
distributions for \ = 0.70 pm, from Deirmendjian (Ref. 5). 

The lower curves represent scattering angles from 0 to 180 
deg, as shown along the bottom of the figure. The upper 
curves depict the small angle part of the same curves on an 
expanded scale, shown at the top of the figure. 
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Figure Captions (Cont. ) 


Fig. 10 


Fig. 11 


Fig. 12 
Fig. 13 


Fig. 14 


Fig. 15 


Fig. 16 
Fig. 17 


ERTS-1 MSS Band 5 image produced from digital tapes by 
the investigators. A 12 x 24 mile area in the western part 
of the Los Angeles basin is shown from the 21 Oct 72 pass. 
The area shown was quite hazy at the time of the ERTS-1 
pass. White lines are included to precisely locate specific 
frame lines in order to correlate the digital analysis with 
the visual image. 

ERTS-1 MSS Band 5 image produced by the investigators 
from digital tapes of the 8 Nov 72 pass. The area shown is 
roughly the same as that in Fig. 10, The area was quite 
clear (haze free) at the time of the ERTS-1 pass. The 
white strips are included to precisely locate specific frame 
lines. 

Similar to Fig. 10 except for MSS Band 7. 

Intensity plots for all MSS bands along line 1871, from the 
same frame of data as Fig. 10. Lines 1865 and 1875 are 
marked on Fig. 10, thus line 1871 can be accurately located. 
The east-west extension, i. e. , point numbers is the same as 
^ or Fig. 10, MSS bands 4, 5, 6 and 7 are plotted in order 
from bottom to top. 

Intensity plots for the line (1788) of the 8 Nov 72 data that 
most nearly coincides with the line plotted in Fig. 13 from 
the 21 Oct 72 data. There is an along scan shift, because 
the frames from these two passes do not exactly overlap. 
These plots are from the same frame of data shown in 
Fig. 11. Lines 1780 and 1790 are marked in Fig. 11, 
therefore line 1788 can be accurately located. 

Histograms of brightness values from ERTS-1 MSS Band 4 
digital tapes. Data from lines 1783 to 1793 and points 1813 
to 2363 were included from the 8 Nov 72 frame. This area 
can be located on Fig. 11. This area includes ocean, beach, 
harbor, vacant land, city, and parkland. The same geo- 
graphical area was included from the 21 Oct 72 frame, i. e. , 
lines 1866 and 1876 and points 1735 to 2285. This area can " 
be located on Fig. 10. 

Similar to Fig. 15, except for MSS Band 5. 

Similar to Fig, 15, except for MSS Band 6. 
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Figure Captions (Cont. ) 


Fig, 18 
Fig. 19 
Fig. 20. 
Fig. 21 
Fig. 22 
Fig. 23 


Fig. 24 


Similar to Fig. 15, except for MSS Band 7. 

Power, spectra from the same data defined in Fig. 15. 

Similar to Fig, 19, except for MSS Band 5, 

Similar to Fig. 19, except for MSS Band 6. 

Similar to Fig. 19, except for MSS Band 7. 

Image of the difference in point by point brightness of 
two frames. The digital data of 21 Oct 72 minus the 
data of 8 Nov 72, covering an area of approximately 
45 x 100 miles, is pictured. Areas of zero difference 
are imaged as middle gray. The most prominent 
features are clouds. Clouds present on 21 Oct' appear 
white, while those present on 8 Nov are black. The 
coastline is visible because the surf was more prominent 
(brighter) on 8 Nov than 21 Oct. The slightly brighter 
area corresponds quite well with the region observed to 
be hazy on 21 Oct. 

Drawing of the essential elements of the solar aureole 
monitor. 


Fig. 25 


Photograph of the solar aureole monitor in operation. 
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Figure 8. Flow chart of the ERTS1 program. 
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APPENDIX A 


GROUND TRUTH DATA 


Table A-I. Visual Observations of the Haze Conditions, 
Including the Altitude of the Top of the Haze Layer Over 
the Test Area at the Times of ERTS-1 Passes 


Letters and numbers in parentheses refer to locations 
shown in Fig. 1. 

In order of increasing haziness the haze conditions are 
described as: clear, very light, light, moderate, 
moderate to heavy, heavy, very heavy, extremely heavy. 

14 Sept 1972 - Pass over San Bernardino. Cloud free. Moderate 
to heavy haze (smog) over the Los Angeles basin (A,B,D, E, H, 

I, L, 1, 6, 7). Haze top at 2, 500 ft. over Torrance (7) and 
Inglewood (H). Heavy haze (smog) over San Bernardino and 
Ontario (4). Clear north of the San Bernardino mountains in 
Antelope Valley and Lucerne Valley. Clear at Big Bear and 
Arrowhead (high in the San Bernardino mountains). 

2 Oct 1972 - Pass over San Bernardino, Broken clouds, 
bottoms around 2, 500 ft, tops about 5, 500 ft. over the Los 
Angeles basin. Very light haze tops at 5, 500 ft. over the 
Los Angeles basin (A, D, E,H, 1,6). Too cloudy over San 
Bernardino and Ontario (4) for analysis. 

21 Oct 1972 - Pass over Los Angeles. Heavy haze (smog) over 
the entire Los Angeles basin (A,B,D, E,G,H,I,L. 1,4, 5,6,7). 
Haze tops at 3, 500 ft. over Torrance (7) and LAX (1), slightly 
higher inland. There is moderate to heavy haze (smog) over the 
southern part of the San Fernando Valley decreasing to light haze 
at the northern side of the Valley. There is a clear area around 
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the Van Norman reservoir. It is clear in the Newhall 
area (J). It is clear over Palmdale (8) and Lancaster (K). 

There are a few small low clouds over the Palos Verdes 
Peninsula (near 7), Orange County (5), and Ontario (4). 

8 Nov 1972 - Pass over Los Angeles. A very clear day. 

It is clear over the Los Angeles basin west of Whittier (I) and 
Arcadia (B). There is very light haze over Orange County 
(5). Light to moderate haze (smog) and scattered clouds over 
Pomona (G) and Ontario (4). It is clear over the San Fernando 
Valley (F, 2, 3). It is clear around Newhall (J) and over 
Palmdale (8) and Lancaster (K}. It is very cloudy northwest 
of Newhall (J). 

25 Nov 1972 - Pass over San Bernardino. There are scattered 
cirrus clouds which might well confuse a haze analysis. There 
is light haze, which tops at 1, 000 ft. over the Los Angeles 
basin. The haze is patchy around Anaheim because it settles 

in low areas. Clear over Riverside and San Bernardino. Light 
haze over Pomona (G) and Ontario (4) similar to the Los Angeles 
basin. 

26 Nov 1972 - Pass over Los Angeles. Smog has been blown out 
to sea. Clear over the San Fernando Valley. Clear over the 
Los Angeles basin except for a very thin layer of smog, 100 to 
200 ft. thick. In some places this layer is at ground level, in 
other places as much as 1, 000 ft. above ground level. It looks 
slightly hazier over Orange County (5) than over the Los Angeles 
basin. 

14 Mar 1973 - Pass over Los Angeles. A very clear day. Just 
a hint of haze over Anaheim and Long Beach (6). Some clouds 
over the San Gabriel mountains. 
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1 Apr 1973 - Pass over Los Angeles. Very light haze (marine 
type not smog) over the Los Angeles basin (A, B, D, E, H, I, 1 , 

6, 7). Haze tops at 5, 500 ft. A few clouds also top at 5, 500 ft. 

No cirrus above. There are heavy clouds northwest of Newhall 

(J). 

19 Apr 1973 - Pass over Los Angeles. Very light haze along the 
coast (D, 1,7) and over the San Fernando Valley (F, 2,3). Moderate 
haze and some low scattered clouds over the remainder of the 
Los Angeles basin (A, B, E, H, I, L, 6). The top of the haze layer 
is at 4, 000 ft. No cirrus. 

7 May 1973 - Pass over Los Angeles. Heavy haze (fog and smog) 
over Los Angeles fca sin (A, D, E, H, 1,7,6), Top of haze layer 
and some scattered clouds at 2, 800 ft. Light fog over the ocean. 
There are cirrus north of Santa Monica (D), but none south of 
Santa Monica. There is moderate to heavy haze over the San 
Fernando Valley (F, 2,3). 

18 July 1973 - Pass over Los Angeles. Moderate to heavy 
haze (fog) over the Los Angeles basin. The haze layer tops 
at 2, 100 ft. Broken clouds top at the same altitude. There 
are high cirrus clouds over the area. 
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Table A -II. Airport Visibility and Humidity Data, Air Pollution 
Control District Data, Vertical Temperature Profile 
and Solar Aureole Brightness 


DATE : 9 ACJ < 3 1972 TIME: 10 : 00 AM P.5.T. 
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Table A-II (Cont. ) 
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Table A-II (Cent.) 
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appendix b 


SOLAR AUREOLE MONITOR 

The solar aureole monitor measured the aureole brightness 
at angular distances of 1.5, 3.0, and 6.0 deg from the center of the solar 
disc along the almucantar, i. e. , along a line parallel to the horizon, 
in a spectral band extending from 0. 6 ^/m to 1. 0 fim. 

A sketch of this instrument is shown in Fig. 24. The incoming 
sunlight passed through the entrance aperture, A. A pinhole irra ge of the 
solar disc passed through the solar disc aperture, O, into the light trap. 
Pinhole images of the aureole were detected at angles of 1.5, 3, and 6 
deg from the center of the solar disc by the detectors, which were 1 mm 
diameter silicon diodes. A W ratten #29 filter was mounted on the face of each 
diode to block wavelengths less than 0. 6 / jm . 

- 4 — (y 

Typical aureole intensities are in the range of 10 to 10 
of the brightness of the solar disc. Thus, the primary technical difficulty 
in constructing an aureole monitor is insuring that a negligible amount of 
the large flux of the solar disc scatters or diffracts into the aureole 
detectors. A pinhole was chosen for the entrance aperture rather than a 
lens, because a lens at that position would need to be inherently a very 
low scatter fabrication and would need to be kept extremely clean. 

Diffraction of light from the pinhole into the aureole detectors can be 
calculated, and thus designed to be negligibly small. The ratio of 
diffracted to incident light is given by (Ref. 12) 



where 


x 


7TD Sin 0 
X 


B- 1 



Jj is the first order Bessel function 

D is the diameter of the aperture 

0 is the angle of diffraction 

A is the wavelength of light 

1/ I o has numerous maxima and minima* Its maximum value falls below 
10‘ 6 for x a* 150, and continues to fall as x increases. 

For 

D = 0. 095" = 2. 413 x 10 3 W m 

X =0.75 ym 

x = 150 

6 =51 min. 

The solar disc subtends a half angle of 1 6 min. This is smeared by the 
finite size of the entrance aperture to 21 min. at the detector plane. 

Thus light diffracted into the 1-1/2 deg detector is much less than 10 ^ 
even from the edge of the solar image, next to that detector. The light 
diffracted into the 3 and 6 deg detectors is much less even than that 
scattered into the 1-1/2 deg detector. Therefore, diffraction from the 
entrance aperture is negligible compared with the aureole signal in this 
system. 

The solar disc aperture, O, allows the intense flux of the 
solar image to pass into the light trap. The diameter of the solar disc at 
this point is 0. 28 inch. This is blurred to a diameter of 0. 375 inch by 
the finite size of the entrance aperture. It is smeared even farther to about 
0. 40 inch by diffraction from the entrance aperture. The diameter of the 
aperture O of 0. 52 inch allows some margin. A separate unit, boresighted 
with the aureole monitor, was used to keep the solar image centered in 
the solar disc aperture. This also maintained the aureole detectors at 
the proper angles. The design of the light trap, baffles, etc. , limits the 
light from the solar disc that scatters inside the instrument into the 
aureole detectors to -2 x 10 of the intensity of the image of the solar disc. 
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A separate unit, boresighted with the aureole monitor, was 
added to measure the brightness of the center of the solar disc in the 
same wavelength band as that used to monitor the aureole brightness. 
Figure 25 shows a solar aureole monitor in operation. 
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APPENDIX C 


IBM BYTE TO CD C WORD CONVERSION PROGRAM. 


This routine for use on Control Data 6000-7000 series computers 
converts information from records generated on 8 bit byte oriented computers 
such as XDS computers or IBM 360-370 computers to convenient CDC 
formats. The routine is written in the CDC assembler language, COMPASS, 
and is compatible with CDC RUN Fortran. 

Depending on the source of bytes for the conversion to CDC 
format the routine may be used in any of three modes. 

Mode 1 Bytes are supplied by the user in the argument array. 

Mode 2 Bytes are supplied by the user in a separate array. 

*Mode 3 Bytes are supplied by the user as blocks on a file. 

The different modes are indicated to the routine by the values in the 
single array which is the argument of the subroutine for all calls. The 
specific values will be discussed later. 

The specific type of conversion desired is indicated by the 
entry point of the routine which is called. These entries are: 

NXBYTE This entry returns bytes to the user as directed by 
the argument array. Bytes are right justified with 
zero fill. The bit pattern is exactly that appearing 
in the desired bytes. No conversions are made. 


'in mode 3 the file reading is done by the routine employing the Record 
Manager under the CDC SCOPE 2.0 operating system on the 7600 computer. 
Use of Mode 3 under different operating systems might require some modi- 
fications to the routine. Modes 1 and 2 are believed to be independent of 
operating system. 
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IHEX 


The CDC display code corresponding to the 
hexidecimal equivalent of the bytes desired is 
returned right justified and blank filled. 

IEBCDC The CDC display code corresponding to the byte is 
returned right justified with blank fill. Upper and 
lower case letters are not distinguished. For the 
assignment of symbols in the CDC character set to 
the IBM hex codes and symbols see Table C-I. The 
symbol v' (display code 66g) is returned for all bytes 
which do not correspond to a letter, a number, or 
one of the symbols defined in Table C-I. 

INT2 The 1 1 s complement integer corresponding to the 16 

bit 2’s complement integer provided in two bytes will 
be returns d. 

INT4 The l's complement integer corresponding to the 32 

bit 2's complement integer provided in four bytes will 
be returned. 

REAL4 The CDC floating point number corresponding to the 
32 bit IBM floating point number contained in four 
bytes will be obtained and returned. 

REAL8 The single precision CDC floating point number obtained 

from the 64 bit IBM format floating point number con- 
tained in 8 bytes by rounding the 56 bit IBM fraction to 
48 bits and converting to CDC format is returned. 

SKPRCD No value is returned. This entry reads a new block 
from the file (Mode 3 operation) and initilizes the 
pointers so that the next call for bytes in sequence 
will start at the beginning of the new block. 

Except for SKPRCD, which returns no values, all entry points 
are used as Fortran functions. All of the entries are called with one 
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argument, an array of at most 10 words. (Mode 1 requires only a 2 
word array. ) The significance of the various words in the array is as 
follows: 

Word 1 The byte location or bytes to be converted. 

If word 2 is positive (mode 2 or 3) -- 

This word indicates the location of the next byte 
desired. 

If word 1 is positive -- 

This word is the number of the first byte desired 
from the array or block for this conversion. The 
first byte of the array or block is numbered 1* 

If word 1 is negative -- 

The next bytes in sequence will be obtained. On the 
first call the first byte of the block or array will 
be returned except in mode 2 operation a negative 
value is not allowed on the first call. 

In both cases this word will be set to -1 on return to 
the user. 

If word 2 is negative (mode 1) -- 

Word 1 contains the bytes to be converted to CDC 
format right justified with zero fill. The contents 
of word 1 will not be changed. 

Word 2 The mode of operation and the number of bytes desired 
for conversion. 

For NXBYTE this must be between 1 and 7, otherwise 
only one byte will be returned. 

For IHEX if this word is positive the bytes are taken from 
the block or array being converted and the value of this 
word should be brtween 1 and 5. If it T s value is 6 or 7, 
then 6 or 7 bytes will be extracted from the array or 
block but only the 6th or 7th byte will be converted. For 
values greater than 7 only one byte will be extracted 
and converted. 
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If the value of this word is between -1 and -5 the 
right hand number of bytes from word 1 of the argument 
array will be converted. If it is zero or less than -5 
only the rightmost byte of word 1 of the argument array 
will be converted. 

For REALS if this word is positive it T s value will be 
ignored and it will be set to 7 on return. 

If this word is negative the most significant byte 
of the 8 required (the IBM sign- exponent byte) will be 
taken as the righthand 8 bits of the complement of this 
word. The other 7 bytes will be taken from the right- 
hand 56 bits of word 1 of the argument array. 

For all other entries only the sign of this word is important. 

If the sign is positive the required bytes are drawn from 
the block or array being converted and word 2 will be set 
to the number of bytes required for that particular 
conversion. 

If the sign of word 2 is negative the value of 
word 2 is ignored and the number of bytes required for 
the conversion are taken from word 1 of the argument array. 


if bytes are being extracted from a block or array (word 2 positive, 
mode 2 or 3 operation) and if bytes beyond the end of the block or array are 
required by the combination of values appearing in words 1 and 2 then mode 
3 operation is assumed and an attempt will be made to read a new block into 
the array. 

The rest of the words in the argument array are needed only for 
mode 2 or 3 operation. See the discussion of the modes to determine which 
should be set by the user in each of the modes of operation. 
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Word 3 
Word 4 
Word 5 


Word 6 


Word 7 
Word 8 


Word 9 


This word contains the length of the array or block 
in bytes* 

This word contains the number of bytes remaining in the 
array or block. 

This word contains the bit position in the word 
corresponding to the next byte to be obtained. This 
quantity is maintained by the routine and should not be 
modified externally. 

This word contains the address of the word in the array 
in which the next byte begins* This quantity is main- 
tained by the routine and should not be modified 
externally. 

This word is not used at present. 

This word indicates the name of the file from which 
blocks are obtained for mode 3 operation. It may be 
an integer, n, from 1 to 99 to indicate the file name 
in the usual Fortran convention of TAPEn. It may 
also contain a file name in left justified zero filled 
display code. In either case the file name must be 
mentioned on the program card. No file card is 
required in the control cards. The file is defined 
by the routine to be RT=U, BT=K and RB=1 and so does not: 
have the usual Fortran file definition. 

This word contains the address of the first word of the 
buffer array which contains the block or array to be 
converted in modes 2 and 3 operation. (This may be 
obtained conveniently with the Fortran function LOCF, 

See the examples. ) 

In mode 2 operation this buffer is filled by the user with 
the array he desires to make conversions from. 

In mode 3 operation this buffer is filled from the file 
defined by word 8, one block at a time. The buffer 
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must be long enough to contain the longest block 
expected* That length will be 8 times the maximum 
number of bytes expected divided by 60 and rounded 
to the next highest integer. 

Word 10 This word contains a flag which indicates the status 
of the file after the most recent call to the routine. 
Possible values are: 

0 No change has occurred in the file or buffer, or a 
new block was read in successfully from the file by 
SKPRCD. 

1 Bytes beyond the end of the current array or block 
were requested. A new record was sa tisfatorily 
read in but no data has been extracted from it. The 
word and byte pointers are set so that the next call 
to the routine will start at the beginning of the 
block or array if word 1 of the argument array is 
negative at that time. 

2 An end of section, end of partition or end of 
information was encountered on the read* The type 
of end of data will be indicated by a message in 

the OUTPUT file. 

3 A parity error was encountered on the last read. 

4 All other record manager detected errors detected 
on the last read are indicated by the value 4. The 
RM error number will be indicated in a error 
message put into the OUTPUT file. 

The various modes of operation, which depend on the source of 
bytes for conversion, are indicated to the routine as described below. 

Mode 1 The bytes to be converted are supplied by the user right 

justified in word 1 of the argument array. This mode is 
indicated by a negative sign in word 2 of the argument 
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array. The value of word 2 is ignored except for 
two cases. For REAL8 conversion word 2 is the 
complement of the most significant byte of the word. 

For IHEX conversion word 2 is the complement of the 
number of bytes to be converted. 

Mode 2 The bytes to be converted are extracted by the routine 

from an array supplied by the user. In this mode the 
following words are set, at least initially by the user. 

Word 1 On the first call this word cannot be negative. A 

specific byte location must be specified. Sub- 
sequent calls may specify the byte desired expli- 
citly or may obtain the bytes in sequence by setting 
(or leaving) this word set to a negative value. 

Word 2 This word must be positive. See the word 
description for limitations on its value. 

Word 3 This must be set by the user to the length, in bytes, 

of the array from which the bytes will be drawn. 

This must be done before the first call to the 
routine and need not be reset unless a new byte 
string of different length is put into the array by 
the user. 

Word 4 This must be set by the user, before the first call 
to the routine, to the length, in bytes, of the array 
from which the bytes will be drawn, i. e. , it is set 
to the same value as word 3. It shouldn't be reset 
by the user thereafter unless a byte string of 
different length is put into the array by the user. 

Word 9 This word must contain the address of the first 

word of the array containing the bytes to be converted. 

Words 5, 6, 7, 8 and 10 are not set by the user. 
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Mode 3 


The bytes are extracted from blocks obtained from a 
file. In this mode the following words are set, at 
least initially by the user. In addition the file name 
must be mentioned in the main program card. 

Word 1 This word may have any value allowed in the word 

description on any call. 

Word 2 This word must be positive. See the word 

discription for limitations on its value. 

Word 4 This word should be set to a negative value on 

the first call of the program to assure reading the 
first block into the buffer. Thereafter this word 
need not concern the user. 

Alternatively the first block may be input with the 
SKPRCD function and this word always ignored. 

Word 8 This must indicate the file name from which the 

records will be obtained. See the word description 
for a discussion of possible forms. 

Word 9 This word contains the address of the first word of 

a storage area used to store blocks as they come 
in from the file. The length required is discussed 
in the word description. 

Words 3, 5, 6, 7, and 10 are not set by the user but may 

contain informiation of interest to him as described 
under the various word descriptions. 

Length: The subprogram is 41 5g words long. 

Other Routine s Required: No other routines unavailable with a normal RUN 

Fortran job are called. Several RUN Fortran Library routines 
are called however in Mode 3 operation. 

Error Mes sages: Error Messages are writtin on the OUTPUT file whenever 

a record manager error occurs or when an end of data is encountered in 
mode 3 operation. 
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Notes and Examples 


1) IEBCDC can conveniently be combained with the routine STRMOV 
to transform long strings of EBCDIC characters into display code character 
strings. For example we assume the mode of the subprogram will be 2. 

(User supplied arrays.) Assume that a user written function routine, IBYTES, 
is used to fill the array IAR with the bit pattern to be converted and that 
the function returns as its value the length in bytes of the array filled. 

Assume that in this particular array there is a string of 108 EBDCIC 
characters starting with byte 325. It is desired to transform that string 
to display code, store it in an array ITST and print it out. 

The following coding will fill the buffer array by means of the 
here undescribed routine IBYTES, set up all required words for mode 2 
use, convert the string of characters desired and print it. 


COMMON TAR (51 2) 

0 1 MENS I ON ICONT nO)tITST(ll)»ITTM<2> 

ICONT (3) =" I BYTES ( I AR 1 ”‘1 

ICONT <4 ) = ICONT (31 
ICONT ( 9 ) = LOCF(IAR) 

ICONT nr = '" 325 “ 

• ITST(ll) = 1H 
00 105 1=1*108 

CALL STRMOV < IERC0C< |CONT)V10*1 * ITST, I) 
10v CONTINUE 

WRITE <6»2> ITST 

— 2~F0RMAT~ (* "ITST =' «*'6Al0*/8Xt5A10P 


The results are: 


ITST 


this isa message. 


^NUMBERS - 0123456789076543210 
r(l=* END OF CHARACTER string. 


SYMBOLS 


>.< ( 
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2) For calls to NXBYTE and IHEX the number of bytes to be 
converted may be, within limits, specified. The results of varying the 
number of bytes even beyond the allowed limits is indicated by the following 
coding, again applied to the character string of the preceeding example. 

This example also illustrates the use of a mode 1 call. This example 
presupposes the initialization coding of example 1. 


ICONT (1) = 325 

"ITST { 2 V -5 
00 135 I = 1*9 
ICONT 

TTST (11= NXRYTE (TCONTV 
ITST (3) = IHEX (ITST) 


TTST (3) = lHt* ' t ■ o i ' 

WRITE (6,5) I TST ( 1 ) * TTST ( 3) ♦ ITST (3 ) 

— . — > . j « n r - \/ A A A 


~ c r r, a m/.t 




135 CONTINUE 

TCOjslT ( 1 ) = 325 

DO 140 1 = 1*6 ' 

ICONT (2) =1 

I0UT = IHEX ( ICONT ) 

WRITE" (6 « 5) TOUT * TOUT 


140 CONTINUE 


The results of executing this are: 



ITST 

ITST 

ITST " 

ITST 

ITST 

ITST" 


00000000000000000343 
00000000000000000310 
"00000000000000144742 
00000000000020144742 
00000000030120152305 
00000016134260343705 
00002264032571 152 302. 
03056636110030040360 
”00000000000000000361 


OOOOOOOOE3 33333333333333330536 

OOOQOOOOC8 33333333333333330343 

OOOOOOC9E2 333333333333034 4 0535 

000040C9E2 3333333337 J303440535 : 

00C140D4C5 33330334373304370340 

E2E2C1C7C5 05350 5350 3340 3420 340 

4005E4D4C2 37330440053704370335 

E2406040E0 05353733413337330633 

OOOOOOOOF1" ‘ 33333333333333330634 


55 5S55 5555555555 0536 £3 

55555555555503430344 C8C9 

"55555555053537330344”' ~ E240C9 

55550535033437330437 E2C140O4 

03400535053503340342 C5 E2E2C1C7 

'55555555555555550437 ' _\ 04. 
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3) These examples illustrates the use of calls for number 
conversion. Assume the array of the previous examples has a sequence of 
14 single IBM precision (32 bit) floating point words starting at byte 
number 121. The following coding will extract these numbers and convert 
them to CDC floating point format. The other numerical conversions may 
be applied in exactly the same way, changing only the entry point name 
used to obtain the various conversions. This example presupposes the 
initialization coding of example 1. 


““ICONT ( 1") 21 

DO 110 I=ltl4 
OUT s RfALA ( TCONT) 

WRITE (6« 3) " 0UT»0UT 

3 FORMAT (* OUT = * ? E20 . 1 0 » 5X » 020 ) ; 

110 CONTINUE 

. ... ... 


The results of executing this are: 


— our 


.toooooooooe* 

. OUT 

ST 

-♦ 1 OOOOOOOOOE* 

OUT 

= 

. 149250Q0G0E* 

our 

= 

-.1492500000E+ 

OUT 

= 

0 » 

OUT 

rr 

-.5000000000E* 

OUT 

- ^ -w — — — 

"* 1562500000E- 

OUT 

x 

- • 1 562500000E- 

OUT 

— 

-.1500000000E+ 

OUT 

x 

Y15000G0000E+ 

OUT 

s 

.5397605347E- 

OUT 

X 

-.5397605347E- 

; OUT' 

" X ~~~ ’ 

•7237005146E+ 

OUT 

3 

-.7237005146E+ 


01 17204000000000000000" 

01 ' 60573777777777777777 

03 17274524000000000000 

03 60503253777777777777 

oooooooooooocooooooo 

00 60603777777777777777 

01 17124000000000000000 

01 60653777777777777777 

02 60540377777777777777 

02 17237400000000000000 

78 13144000000000000000 

78 64633777777777777777 

76 231 47777777700000000 

76 54630000000077777777 
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Table C-l. IEBC'DC Svmbol Conversion Code 


HEX 

Value 

IBM Symbol 

CDC Display Code 
(Octal) 

CDC 

Symbol 

40 

space 

55 

blank 

bA 

4 

75 


bB 

• 

57 


be 

< 

72 

< 

bD 

( 

51 

( 

bE 

+ 

b5 

+ 

bF 

1 

71 

1 

50 

& 

70 

1 

5A 

I , 

• 

66 

V 

53 

$ ' 

53 

4 . 

4> 

5C 

* 

^7 

-x- 

53 

) 

52 

) 

5E 

> 

77 

j 

5P 

—\ 

74 

<> 

6o 

- 

he 

- 

6l 

/ 

50 

/ 

Cb 

} 

56 

> 

6 c 


76* 

— 1* 

6d 

— 

65 

n 

6e 

> 

73 

> 

6f 

7 

67 

A 

7A 

Z 

63 

: 

7B 

# 

60 

= 

7C 

@ 

61 

c 

7D 

« 

62 

3 

7E 

- 

5 b 


7F 

n 

6b 

* 

* 

Octal 76 prints % on 200 

terminals and — | on high 

speed Line prin 
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